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A Sawtooth Permanent Magnetic Lattice for Ultracold Atoms 
and BECs 
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Abstract. We propose a new permanent magnetic lattice for creating periodic arrays of loffe-Pritchard 
permanent magnetic microtraps for holding and controlling ultracold atoms and Bose-Einstein condensates 
(BECs). Lattice can be designed on thin layer of magnetic films such as r66GdiO-Fe8oCo4. In details, we 
investigate single layer and two crossed layers of sawtooth magnetic patterns with thicknesses of 50 and 
500nm respectively with a periodicity of Ifim. Trap depth and frequencies can be changed via an applied 
bias field to handle tunneling rates between lattice sites. We present analytical expressions and using 
numerical calculations show that this lattice has non-zero potential minima to avoid majorana spin fiips. 
One advantage of this lattice over previous ones is that it is easier to manufacture. 

PACS. 67.85.Hj ,67.85.Jk,37.10.Gh,37.10.De,85.85.-f j 



1 Introduction 

Focused laser beams have been used extensively in years 
for cooling atoms and molecules to hold them [T1[UI31 
mis]. These optical traps are made by a radiation field 
whose intensity has a maximum in space. Due to inter- 
action with electric field of laser, small clouds of atoms 
can be confined, manipulated and controlled in optical 
lattices [51[71|5]. Radiation absorption and induced dipole 
forces constrain atoms to optical lattice |9j. Periodic op- 
tical lattices are the result of interference of intersect- 
ing laser beams and have allowed for the observation of 
the Mott-insulator to superfluid transition 10,11^, stud- 
ies of low-dimensional and quantum degenerate gases [l2l 
113) and coherent molecules [141 . In quantum computation, 
optical lattices provide a large qubit system and produce 
gate operations on the register [T5l[l6] . 
Using magnetic field is another approach to realize har- 
monic potential in space [T7]. Magnetic lattices are cre- 
ated using periodically magnetized films or current car- 
rying wires on atom chips 18,19,20,21,22 . These mag- 
netic potentials are promising alternatives to the optical 
lattices and provide a high degree of design flexibility, ar- 
bitrary geometries and lattice spacing. Magnetic lattices 
can produce highly stable periodic potential wells with 
low noise j23j. Single 3D magnetic microtraps also have 
been introduced either by current carrying wires or per- 
manent magnetic slabs [24ll25| . Recently, development of 
technology has let us for working on high-quality mag- 
netic materials based on permanent, perpendicularly mag- 
netized films [26u27j . Periodic grooved magnetic mirrors 
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have been made by the technology of thin films to re- 
flect ultracold atoms ^B]. In the presence of magnetic 
bias fields, mirrors are transferred into ID arrays of 2D 
waveguides and 2D arrays of 3D traps ^8] . Magnetic films 
such as TbeGdiOFesoCoi have excellent magnetic prop- 
erties; homogeneity, high magnetization, large coercivity 
(iJc -3kOo) and high Curie temperature(-300°C) [H 
[II. 

In this paper, we introduce a new atom chip based on per- 
manent periodic grooved magnetic film and external mag- 
netic fields. This sawtooth magnetic lattice is proposed 
for controlling, manipulating and trapping clouds of ul- 
tracold atoms when they are in low-field seeking states. 
Analytical expressions and numerical calculations show 
that the magnetic field minima are non-zero. Therefore, 
ultracold atoms can be hold and trapped without any spin 
flips in lattice sites. In particular, we investigate ID array 
of this grooved magnetic lattice, when bias fields are ap- 
plied. Results show that produced waveguides have con- 
trollable depth and frequencies to handle clouds of quan- 
tum degenerate gases on the chip surface. We show that 
two crossed separated layers of periodic arrays also can 
produce 2D magnetic lattice of 3D microtraps when bias 
fields are applied. The microtraps have non-zero minimum 
and high depth (~0.2mK) which cause prevention of ma- 
jorana spin flips and efficient restriction of trapped atoms, 
respectively. Magnetic atom chips based on permanent 
magnetic materials have some advantages over optical or 
microwire devices. Since there is no spontaneous emission 
due to decay of excited atoms, light scattering and deco- 
herence of cold atoms vanishes and no beam alignment 
is required [5]. Heat dissipation and loss of atoms due 
to fluctuation of current in current-carrying microwires 
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are prevented in permanent magnetic lattices. They can 
produce large field curvatures, large trap depth and high 
frequencies with very low technical noises [23^. the Main 
motivation of this work is introducing a more flexible per- 
manent magnetic lattice compared with the configuration 
of the two crossed arrays of rectangular permanent mag- 
netic slabs proposed in [18 . Considering the shape of this 
new atom chip, current technology in grooving of films 
for magnetic mirrors allows us to produce it more easily 
compared with previous ones. 



2 Analytical expressions for Magnetic lattices 

We study first the simple case of a single layer of sawtooth 
periodic array of grooved magnetic films with a thickness 
of t, periodicity a along the y-direction, magnetization 
M=Mz'z and constant bias magnetic fields Bix, Biy and 
Biz along the x-, y- and z-directions. 



2.1 Atom waveguide 

The configuration of the permanent magnets is shown in 
Fig.l. Pattern of magnetic film has a periodicity (distance 
between two alternative peaks) of a = 1/im and a thick- 
ness oi t = 500nm while dimensions of lattice in x- and 
y- directions, Ix and ly are 1mm (Fig. 1(c)). For analytical 
and numerical calculations, we approximate our sawtooth 
magnetic pattern thin rectangular magnetic slabs same 
magnetization and periodicity but with different widths. 
As Fig.l (a) shows, infinite periodic arrays of magnets are 
placed above each other to shape sawtooth structure. At 
a large distance where k{z — t)^l, the using expressions 
given in Ref. |18j . the components of the magnetic field 
can be written 

= Si, (1) 
By = Bq^ sin (fcy)e"''"' + Biy 
B^ = Bo^ cos {ky)e-^'- + Bi^ 

where k=2Tr/a and Bq^ — Bq'^^^^ ji, respectively. 
Bq and 7i are Mz/7r(Gaussian units) and related coeffi- 
cient of the ith infinite periodic magnets on atom chip 
[see Fig. 1(b)], respectively. 7i is determined from the se- 
ries [for more details, see section 2.3] 

^ {2j + l)ka 

where n is the number of sub-layers (i.e. order of ap- 
proximation) . The magnitude of the magnetic field above 
surface can be written 



(a) 




± 



Fig. 1. (Color online) (a) ID array of sawtooth permanent 
magnetic lattice, (b) Schematic of sawtooth permanent mag- 
netic lattice film which has been estimated by n=10 sub-layers 
of rectangular magnetic slabs. Periodicity of the sawtooth per- 
manent magnetic lattice is the periodicity of the rectangular 
slabs and is defined by distance of two alternative peaks of the 
structure. Lattice has the periodicity a—lfim, magnetization 
Mz=3.8 kG along the z-direction and thickness t=500nm. For 
n—10, thickness of each sublayer is t/n—50nm while width of 
slabs is changed according to Wi= a{i+l)/n. (c) Ix and ly, di- 
mensions of the atom chip in x- and y- directions, respectively. 
For nr.=1000 sawtooth magnetic slabs 1^= Zj; = lmm. External 
magnetic field is applied along x- and {/-directions to produce 
waveguides with non-zero potential minima. 
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(a) 



B(y,z) = {Bf^ + Bfy + Bt,^ 
B cos {ky)]e^''^ 



2[Bo^BiySin(fcy)+ (3) 



In simple case when we take Biz=0, the minimuni of 
magnetic field in ID array of periodic magnetic lattice of 
2D microtraps is 

Bmin = I Six I (4) 



where the center of the waveguides are located at 



(5) 



\B 



lyl 



where Uy is the trap number in the y-direction and takes 
values of 0, ± 1, ± 2,... . 

According to equation 4, the x component of the bias field 
is necessary to avoid spin flips near the center of 2D mi- 
crotraps. The magnetic barrier heights along y- and z- 
directions which depend on the bias fields are given by 



ABy 

AB'- 



max 



Bn 
Br. 



{Bl^+AB%f''~\B,^\ (6) 
{Bl^ + Blyf/^^\B,^\ 



Second derivatives (curvatures) of the magnetic field at 
the minimum of the potential and trap frequencies of the 
waveguides along y- and z-directions could be analytically 
determined by 



d'^B d^B 



dz 



r2 



Six I 

mSix 



^l/2 



(7) 
(8) 



where /i^, g^, m-^ and m are the Bohr magneton, 
Lande factor, magnetic quantum number and mass of atom, 
respectively. Equations (5) and (8) show that increasing 
Biy reduces z 




Fig. 2. (Color online) (a) Two layers of sawtooth configuration 
permanent magnetic lattice . (b) Schematic of (a) estimated by 
n=10 sub- layers of rectangular magnetic slabs for each layer. 
Top and bottom arrays have corresponding coefficients series 
\i and 7i, respectively in magnetic field. Periodicity of top 
and bottom arrays are a=l^xa along both x- and y-directions. 
Arrays have magnetization Mz=3.8 kG along the z-direction 
and thickness fi=500nm and t2=50nm for bottom and top 
arrays, respectively. 



components of the magnetic field, for the large distances 
(k(z-ii-<2)^l) from top of the surface can be obtained by 

B^ = Box sin {kx)e-''' + Bi^ (9) 

By = Bq^ sin {ky)e^''^ + Biy 

Bz ^ (B07 cos {ky) + Box cos {kx))e-''^ + Bi^ 

where Sqa ~ BqJ^^i- Coefficient of the ith sublayer of 
infinite magnets of the top array, is determined, simi- 
lary 



xions ana [S) snow mat mcreasmg ^ ' (2j + l)fca, , ^t, k((i-^)t +t 

and increases frequencies of the traps, Ai = ^2^*^°^' ■''^ " ~ ^ > ^ 



respectively. The frequencies of traps are independent of 
magnetization of the magnetic film. They depend on the 
external magnetic fields and the periodicity of magnetic 
lattice. 



2.2 2D permanent magnetic lattice of microtraps 

In order to produce loffc-Pritchard microtraps, two crossed 
layers of periodic arrays of sawtooth magnetic films with 
two different thicknesses of ii and ^2, separation of d and 
same magnetization Mz with uniform bias fields is inves- 
tigated (Fig. 2). The top and bottom arrays have same 
periodicity a along x- and y-directions. Thickness of bot- 
tom array, ti is larger than thickness of top array, ^2- 
is along z-direction and bias fields are Bix and Biy. The 



+d) 



(10) 



where d is the gap between bottom and top arrays. 
The magnitude of magnetic field above the surface can be 
written 

B(x, 7/, z) = {Bl^ + Biy + Bl + 2[BoaSix sin {kx) 
+BQjBiy sin {ky) + BqxB^ cos (kx) 
+Bn^B, cos (fcy)]e-'=" + [B^^ + 
+2BoxBa-y cos (kx) cos (fc2/)]e~2fe^)i/2 (i^) 

Non-zero minimum of the magnetic field and the center 
of microtraps in x-,y- and z-directions for Biz=0 can be 
determined by 
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Vn 



\Bi^Bo\ — i3iy_Bo^| 

In 

(n^ + -)a 



(12) 
(13) 



where rix and riy arc number of microtraps and take 0, 
±1, ±2,... values. Diniensionless constant a is (1/2)(1+1/q;o) 
where ao — Box/Bo-y- At the center of 3D microtraps, the 
curvature of magnetic fields and barrier heights can be 
obtained by 



d^B 


d^B 


Id^B 






2 dz^ 




= ABy = 


a2\Bi^\ 


AB'- 


= (1 + «o 


Y"\B.A 



= fc'ai|Bix| 



(14) 
(15) 



where ai and a2 are 



ai 



2al 



(l + a2)i/2|l-«o| 
The frequencies of microtraps in three directions are 



(16) 



^ ^//fB5F^;^V^ixai .1/2 
y/2 m ' 



(17) 



For a symmetrical 2D magnetic lattice, same magnetic 
height barriers in x- and y-directions, AB'-^ = ABy are 
needed. The constraint of Biy ~ a^Bix satisfies our re- 
quest where uq — Bqx/Bqj. 



2.3 Corresponding magnetic coefficients of estimated 
layers 

The components of magnetic field for rir infinitely long 
and thin permanent magnetic slabs along the a;-direction 
arranged with periodicity a along the y- direction and sep- 
aration distance a/2 between the center of neighboring 
slabs, are given by [29] 



(18) 



Bx = 



Boy sin {ky)e 
-Bz = Bay cos {ky)e 



-kz 



-kz 



where k and Bfjy are 27r/a and i3o(e'^*-l), respectively. 
As we mentioned earlier in 2.1 and 2.2, permanent mag- 
netic slabs with sawtooth geometry can be estimated by n 



a/2 

Fig. 3. (Color online) Schematic of one slab in permanent 
magnetic array with periodicity a and width a/2. Slabs divided 
into arbitrary (e.g. n=10) narrow slabs A with same period- 
icity. The magnetic field of narrow slabs related to initial one 
by coeflicient ci. In Eq.l9 to obtain ci, slab A shifted to right 
and left sides by small steps z4(i-f 1/2) , where i=0, 1,..., n/2- 
1. Summations are taken on all dashed spaces and results are 
compared with magnetic field of large slab to find Ci. 



multi-layer arrays of rectangular magnets with same peri- 
odicity a and different width Wi , where i is the number of a 
specific sublayer. The major problem is the calculation of 
magnetic field of sub-layers with Wi >a/2 or Wi <a/2. In 
other words, for magnetic slabs with different width and 
separation, Eq.l8 does not work and needs to be modified. 
As Fig.l (b) shows, width plus separation between neigh- 
boring slabs for each sublayer i is constant and equal to 
periodicity of the lattice. It means that for any modifica- 
tion of Eq.lS, k remains unchanged. The new components 
of magnetic field for slabs with arbitrary width are related 
to Eq.18 by coefficient 7. The slabs of periodic lattice by 
width and separated gap a/2 can be divided into n nar- 
row slabs (Fig. 3) with width A=a/2n. The magnitude of 
magnetic field due to these narrow slabs is related to the 
initial slab with coefficient ci. For the calculating of ci, 
narrow slabs can be shifted by small steps of a /An along 
y axis to make initial array. The component of the mag- 
netic field for these shifted narrow slabs, for example in 
y-direction must be consistent with Eq.lS, so we have 



n/2-1 

Ba^v = ciBoy J2 sinfc(y ± (2j + l)A)t 

3=0 

= Boy suv{ky)e~^'' 



(19) 



where ci is determined 



Cl 



n/2-] 



cosfc(j/ + (2j + l)kA)y 



(20) 



We use these narrow slabs A as unit to design any periodic 
lattice by arbitrary slab width Wi and periodicity a. In 
general, the A slabs must be several times shifted to the 
right and left to model the given periodic array. These 
shifts lead to another series in calculation of the magnetic 
field which is given by 



7, = - V cos(2j + l)kA{e^^ - l)e^'^^-'^) (21) 

3=0 

For sawtooth array, all estimated sub-layers have same 
thickness t/n and distance from origin l-(i-|-l)/n. 
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Summation on 7^ and for n=10 sub-layers are deter- 
mined 13.97 and 5.47, while for n—20 sub-layers reduce to 
13.02 and 5.45 values, respectively. 



3 Numerical calculations for magnetic lattice 

Numerical calculations for finite arrays of sawtooth mag- 
netic films with finite length were made. Radia Package [3U] , 
interfaced to Mathematica software was used to numeri- 
cally obtain magnetic field for different configurations of 
magnets in ID array of waveguides and 2D arrays of 3D 
microtraps. 



3.1 Atom waveguide 

Numerical calculations of ID array of waveguides above 
the sawtooth surface of the atom chip were performed 
when bias magnetic field along the y-direction, Biy was 
applied. Thickness t =500nm, periodicity a=l/xm, mag- 
netization 47rAf2=3.8kG, finite toothy magnets number 
71^=5001 which is approximated by n—lQ rectangular slabs. 
In Fig. 4, the components of magnetic field of atom chip 
above the surface are shown where Biy= — 5G. This con- 
figuration may be used as spatial grating to diffract slowly 
moving atoms. By including Bix= — 3G parallel to the 
waveguide axis, center of waveguides becomes non-zero 
with Bjnin=\Bia:\=SG to avoid spin flips of cold atoms. 
Center of waveguides is at 2;mi„=1.29/im from bottom 
of the layer. Magnetic barrier heights along y- and z- 
directions are 7.4 and 2.8G, respectively. The trap frequen- 
cies of 2D microtrap in y- and z-directions are numerically 
determined 27rx36.7 kHz for ^''Rb atoms. Fig. 4 (a) and 
(b) show excellent agreement with numerical calculations 
and analytical expressions of magnetic field along y- and 
z-directions for z — t 3>a/27r. Numerically calculated pa- 
rameters and related analytical expressions for a sawtooth 
magnetic film with bias fields are listed in table 1. 
Large curvatures of the magnetic field (~10^° G/cm^) can 
be realized in the permanent magnetic lattices due to their 
high magnetization which a small size of the microtrap is 
led to high frequencies for efficient trapping and control- 
ling of ultracold atoms. Rate of spin flips transition for 
single atom in ground state of harmonic waveguide, when 
nuclear spin is neglected and total spin is spin of single 
electron, can be written [33] 



r = 



(22) 



where Umin and uj are potential energy minimum and 
radial oscillation frequency, respectively. Eq.22 works for 
Umin ^ ff-^- By using table 1, where J7mi„/fcB=201.5/ik 
and haj/kB=l-7^k, decay life time of trapped ^^Rb atoms 
in F=2 and mp=2 state are so long and there are no ma- 
jorana losses in the waveguides. Thereby, the introduced 
waveguides can be used to transfer cold atoms over the 
surface of single layer atom chip along x-direction with- 
out any significat atom losses from the harmonic trap. 



12 
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Fig. 4. (Color online) (a) and (b) The magnitude of mag- 
netic field close to the center of the waveguides in y-and z- 
directions, respectively. Numerical calculations (red solid) and 
analytical expressions (dashed-blue) are obtained using Radia 
and Eq.3, respectively, for n—10 estimated sub-layers. Plots 
show good versatility between numerical and analytical deter- 
minations for n, =5001 sawtooth slabs with magnetic bias field 
Bmin~\Bix\= 3G and Biy=—5G. Center of the waveguides is 
at Zmin = l-29iim and ymin = {ny + l/4)a, 
waveguides. 



3.2 2D permanent magnetic lattice of microtraps 

Here, we investigate two crossed sawtooth magnetic lay- 
ers to create 3D lofFe-Pritchard microtraps. This magnetic 
lattice is created when bias fields are applied. Fig. 2 shows 
structure of the atom chip. In each layer, we have consid- 
ered 71^=5001 parallel long magnets of triangular cross 
section with magnetization 47rM2=3.8kG and periodic- 
ity a=l/xm. The bottom and top layers have thicknesses 
ti=500nm and i2=50nm, respectively. The number of sub- 
layers in each layer is n =10 and there is no separation 
between them {d—0 in Eq.lO). Each sawtooth magnetic 
layer can be estimated by finite number of sublayers of 
parallel rectangular magnetic slabs with same periodicity, 
width Wi={i+l)a/n and distance Si—{l-{i+l)/n)t from 
the plane z — 0, where i—0, 1,..., n-1 is the number of 
the ith sublayer. To create crossed permanent magnetic 
sawtooth layers, according to Fig. 2(b), Si is subject to 
conditions 0< Si < ti and ti < Si < t2 for bottom and 
top arrays, respectively. 

For n=10 sublayers, great versatility can be seen in Fig. 5 
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Table 1. . Numerical and analytical determinations for single array of sawtooth permanent magnetic lattice. Lattice has 
periodicity o=l/um along y-direction, thickness /;=500nm and magnetization Mz=S.8 kG along ^-direction (Fig.l(b)). Results 
is obtained for ^^Rb atoms in different approximation orders (sub-layers) with bias field Bix=—SG and Biy=—5G. 

Parameter Definition Numerical Analytical 

n=W n=20 n=W n=20 

rir Number of magnets in j/-direction 5001 oo 

yminifi^) Center of first waveguide in j/-direction 0.25 0.25 





Center of waveguides in direction 


1.29 1.28 


1.29 1.28 


#f (^) 


Curvature of B in j/-dircction 


3.29x10^° 


3.28x10^° 


d^B / G \ 


Curvature of B in z-dircction 


3.29x10^° 


3.28x10^° 


c^y/27r(kHz) 


Trap frequency in y- direction 


231.5 


231.5 


w^/27r(kHz) 


Trap frequency in z- direction 


231.5 


231.5 


Bmin (C) 


Minimum of B at center of waveguides 


3 


3 



Table 2. Numerical and analytical determinations for two crossed layers of sawtooth permanent magnetic lattice. Lattice 
has periodicity a=l/im along x- and y-directions, with bottom and top thicknesses ti=500nm and t2=50nm respectively. Two 
layers have same magnetization M2r=3.8 kG along 2;-direction (Fig.2(b)). Results are obtained for ^'^Rb atoms and bias field 
Bix=—5.6G and Biy= — 2.18G for n=10 sublayer and Six=— 5.6G and Bij,=— 2.3G for n=20 sublayer. 



Parameter Definition 



Numerical 
n=10 n=20 



Analytical 
n=10 n=20 



rir 

Xmin {lJ.ni) 
J/min(/im) 

(^) 

(^) 

W:„/27r(kHz) 
Wj;/27r(kHz) 
W;^/27r(kHz) 

Bmin (C) 



a^B / G 
d^B / G 



Number of magnets in x- and y-dircction 5001 

Center of first microtrap in a;-dircction 0.25 

Center of first microtrap in j/-direction 0.25 

Center of microtrap in ^-direction 1.33 1.31 

Curvature of B in a-direction 2.98x10^° 

Curvature of B in y-direction 2.98x10^° 

Curvature of B in z-dircction 5.95x10^° 5.97x10^° 

Trap frequency in x- direction 220 220.5 

Trap frequency in y- direction 220 220.5 

Trap frequency in z- direction 311.6 312 

Minimum of B at center of microtraps 4.41 4.27 



oo 
0.25 
0.25 
1.33 1.31 
2.98x10^° 
2.98x10^° 
5.96x10^° 5.97x10^° 



220.3 
220.3 
311.6 
4.41 



220.5 
220.5 
312 
4.27 



(a), (b) and (c) between analytical and numerical plots of 
magnetic field along the x-, y- and ^-directions, close to 
the center of the permanent magnetic lattice. By adding 
more sublayers, we expect to have a sawtooth structure 
with better accuracy. Numerical and analytical results for 
the 2D magnetic lattice have been compared in tables 2 
and 3 for n= 10 and 20, respectively. For example, sym- 
metrical magnetic traps in x- and y- directions for n=10 
can be obtained at Zmin=l-33y^m from bottom of the lay- 
ers by applying the bias field Bix= — 5.6G and Biy= 
— 2.18G. The minimum of the magnetic field at the center 
of each microtrap is 4.41 G. 
According to Fig. 6, a bias magnetic field along the x- 
direction can change the trap frequencies. It is also pos- 
sible to create permanent magnetic lattices with different 
barrier heights in the x- and y- directions. For example, 
if we apply a bias magnetic field with Bix= — lOG and 
Bij,=— 0.87G the trap depths along the x- and y-directions 
4.5G and 15G, respectively. For this bias field, correspond- 



Table 3. Numerical results, using Radia, and analytical ones 
for *^Rb in F=2 and mF=2 state. Results have been calculated 
for the structure shown in Fig. 2. 



l-*a.ra.iiictcrs (//k) 


Dt'tiiiititni 


;; = iO 


;;=20 


Umin I^B 


Mininmin Potential 


296 


287 




Potential barrier height 


201 


213 




in x-direction 






AV/kB 


Potential barrier height 


201 


213 




in y-direction 






AWjkB 


Potential barrier height 


107 


119 




in ^-direction 






f)Wx/kB 


Energy level spacing in 


10 


10 




a;-direction 






hWy/kB 


Energy level spacing in 


10 


10 




y-direction 






hjJz/kB 


Energy level spacing in 


15 


15 



0-direction 
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Fig. 5. (Color online)(a), (b) and (c) The magnitude of magnetic field close to the center of the microtraps in x-, y- and 
z-directions, respectively. Numerical calculations (red solid) and analytical expressions (dashed-blue) obtained by Radia and 
Eq.ll respectively, for n=10 estimated layers . Plots show excellent versatility between numerical and analytical determinations 
for nr-=5001 sawtooth slabs with magnetic bias field _Bii: = ~5.6G, -Bij,=— 2.18G and -Bmin=4. 41G. Center of microtraps along 
^-direction is at 2mi„=1.33^m. Position of first microtrap along x-and {/-direction is Xmin=ymin=0. 25 jj,m which is repeated 
by periodicity of arrays, (d) and(e) Contour plots of magnetic field at z—Zmin and x—Xmin planes, respectively. Yellow points 



(grey, in black and white) show positions of microtraps . The density of lattice sites is ^ 
field on the surface of atom chip at z—Zmin plane. 



^10 traps/mm . (f) 3D plots of magnetic 



ing frequencies along the x- and y-directions are not the 
same and atoms feel two different barrier heights in the 
X- and y-directions, respectively. Decreasing or increasing 
potential barrier heights will change the rate of tunneling 
between lattice sites. Different rates is useful in quantum 
tunnelling experiments. 



4 Conclusion 

We have introduced a new permanent magnetic lattice 
based on perpendicularly magnetized thin sawtooth films. 
This lattice can produce loffe-Pritchard 2D and 3D mi- 
crotraps in ID and 2D lattices, respectively for holding 
ultracold atoms and BECs when a bias field is applied. 
We have presented analytical expressions for parameters 
such as location of the microtraps, curvatures of the mag- 
netic field and also frequencies and trap depth close to 
the center of the lattice. We also have shown that nu- 
merical calculations are in agreement with analytical ones 
for two different number of the sublayers n. The compo- 
nents of homogenous external magnetic field along x- and 
j/-directions, Bi^ and -Bit,, can change the trap depth and 




3 3.5 4 4.5 5 

Bx ( G) 

Fig. 6. (Color online ) Plots show effect of the x component 
of the bias field on the trap frequencies for a 2D permanent 
magnetic lattice. 



frequencies of microtraps. Control over the trap depths 
let us to handle tunneling parameters via increasing or de- 
creasing potential barrier heights for cold atoms and BECs 
which has wide applications in quantum information pro- 
cessing and superfluid to Mott insulator quantum phase 
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transition studies. Trap density is ~ 10 atoms/mm for 
loading cold atoms and BECs. Considering the current 
technology, it should be easy to manufacture the proposed 
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